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Abstract 

A representation of the static MHD equilibrium system in coordinates connected 
with magnetic surfaces is suggested. It is used for producing families of non-trivial 3D 
ly-j ■ exact solutions of isotropic and anisotropic plasma equilibria in different geometries, 
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finding coordinates in which exact equilibria can be constructed are discussed; examples 
and their applications as physical models are presented. 
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1 Introduction 



X 

■ In this paper we consider two systems of partial differential equations that are most fre- 

quently employed for continuum description of plasmas. The isotropic Magnetohydrody- 
namics (MHD) system has the form £Q 

|^ + divpV = 0, (1) 

<9V 1 V 2 

p— — = pV x curl V B x curlB — gradP — pgrad h /iiAV , (2) 

at fi 2 

<9B 1 

— = curl(V x B) + r/AB n = — , (3) 

at a\i 

divB = 0, J = -curlB. (4) 
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Here V is plasma velocity, B, magnetic field, J, electric current density, p, plasma density, 
P, pressure, p, the magnetic permeability of free space, a, conductivity coefficient, pi, the 
plasma viscosity coefficient; p, resistivity coefficient. 

For a vanishing magnetic field, B=0, the above system is reduced to Navier-Stokes equa- 
tions of motion of a viscous compressible fluid. 

The MHD system correctly describes the medium maintained isotropic by frequent colli- 
sions. However, when the mean free path for particle collisions is long compared to Larmor 
radius (e.g. in strongly magnetized or rarified plasmas), the Chew-Golberger-Low (CGL) 
model [2] is relevant. Like the MHD system, the CGL equations are derived from Boltzmann 
and Maxwell equations, but the density function in Boltzmann equation is expanded in the 
powers of the Larmor radius. In the CGL model, the gas pressure P is replaced by a 3 x 3 
pressure tensor 

Vij=pJi j + TB i B j , r = ^^, i,j = 1,2,3 (5) 

with two independent components: the pressure along the magnetic field p\\ and in the 
transverse direction pj_. 

Under an assumption p\ = r] = 0, i.e. in the case of non-viscous infinitely-conducting 
plasmas, both MHD and CGL systems have several remarkable analytical properties. In 
particular, one can name, for both systems, the "frozen- in magnetic field" property (Kelvin's 
theorem), Lagrangian and Hamiltonian formulation [3], and conservation of helicity [3]. This 
approximation is natural in the case of large kinetic and magnetic Reynolds numbers, and 
is used in this work. 

The present paper is devoted to the study of equilibrium plasma configurations modeled 
by MHD and CGL equations. The MHD equilibrium system is 

1 V 2 
pV x curl V B x curl B — grad P — p grad -— = , (6) 



div(pV) = , curl(V x B) = , div B = . (7) 

The equilibrium CGL system can be put to the form 

(\ \ V 2 
pV x curl V — t B x curl B = grad p± + p grad h 

V/* / B2 2 (8) 

+t grad — + B(B ■ grad r) , 



div(pV) = , curl(V x B) = , div B = . (9) 

The above systems are closed by appropriately chosen equations of state (one for MHD 
equilibria, two for CGL equilibria.) In this paper we restrict attention to incompressible 
plasmas 

divV = 0. (10) 
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Incompressibility approximation is commonly used for subsonic plasma flows with low Mach 
numbers M< 1, M 2 = Y 2 /(>yP/p). (Here 7 is the adiabatic exponent.) For incompressible 
plasmas, the continuity equation divpV = implies V-gradp = 0, hence density is constant 
on streamlines. 

Both of the systems under consideration are non-linear systems of partial differential equa- 
tions essentially depending on three spatial variables. No general methods of construction of 
exact solutions to the corresponding boundary value problems are available; the question of 
stability is answered only for particular types of instabilities (for a review, see [H].) However, 
some progress have been recently achieved in constructing particular exact solutions, which 
are also demanded by applications. In this connection, one can mention exact solutions 
obtained using reductions by symmetry groups (e.g. Grad-Shafranov and JFKO equations, 
see j7J El HE [Jill UH E2]) and solutions constructed using Euler potentials (e.g. [T3l IT1] ). 

The main goal of this paper is to present a more general method of construction of exact 
solutions to MHD and CGL plasma equilibrium equations and their static reductions, in 
different geometries, and with different physical properties. 

The method of construction of exact solutions described below is based on the intrinsic 
property common to both dynamic MHD and CGL equilibria and many static cases - the 
existence of magnetic surfaces |Hl EH] . 

In many important cases, an orthogonal coordinate system can be constructed, with one of 
the coordinates constant on magnetic surfaces of some plasma equilibrium configuration. In 
such coordinates, the static plasma equilibrium system is reduced to two partial differential 
equations for two unknown functions. One of the equations of the system is a "truncated" 
Laplace equation, and the second has an energy-connected interpretation (Section |21) 

The suggested representation of the static plasma equilibrium system is used for producing 
particular exact solutions for static and dynamic, isotropic and anisotropic plasma equilibria 
in different geometries. 

Formulas presented in Section |3] give rise to explicit expressions for static force- free plasma 
equilibria in coordinates with particular relations between metric coefficients. Also, in many 
classical and non-classical systems of coordinates, non-trivial gradient vector fields can be 
built, tangent to prescribed sets of magnetic surfaces fSection [3.2l ) Though gradient fields by 
themselves represent only degenerate plasma equilibria with constant pressure and no electric 
currents, they can serve as initial solutions in infinite-parameter transformations (such as 
Bogoyavlenskij symmetries [HIE] of Plasma Equilibrium equations, and transformations from 
MHD to CGL equilibria OEICSI, which produce non-trivial equilibrium configurations in 
MHD and CGL framework, with non-vanishing plasma parameters. 

In Section |3J we use the new representation to construct families of MHD and CGL plasma 
equilibria. We start from static gradient and non-gradient solutions with magnetic surfaces 
being nested spheres, ellipsoids, non-circular cylinders, and surfaces of other types. These 
solutions, by virtue of Bogoyavlenskij symmetries and MHD-to-CGL equilibrium transfor- 
mations, give rise to infinite families of more complicated dynamic and static, isotropic 
and anisotropic equilibrium configurations. In the majority of constructed equilibria, the 
behaviour of magnetic energy and other plasma parameters within the plasma domain is 
physical. 
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Auxiliary statements of the Section 13.31 allow, in many cases, the extension of static equi- 
librium magnetic fields with a Killing component, thus changing the equilibrium topology. 
Hence the application of these transformations before Bogoyavlenskij symmetries results in 
the change of the domain of arbitrary functions of the latter. 

The value of some of the suggested solutions as models of astrophysical phenomena is 
discussed. It is shown that some essential features of the models and the relations between 
macroscopic parameters are in the agreement with astrophysical observations. Unlike the 
majority of existing models, the presented solutions are exact and generally non-symmetric. 



2 The representation of Plasma Equilibrium equations 
in coordinates connected with magnetic surfaces. 

The general MHD equilibrium system ©-0; its static reduction 

curl BxB=/i grad P , divB = 0, (11) 
and the force-free plasma equilibrium system (P = const) 

curlB = a(r)B, divB = (12) 

are known to possess a family of 2-dimensional magnetic surfaces (or a foliation) ^(r) = 
const, to which both velocity V and magnetic field B are tangent, and thus magnetic field 
lines and plasma streamlines lie on these surfaces 1 [HI EI] • 

Noting the value of representation of the static Plasma Equilibrium system ()11|) in special 
coordinates (e.g. the derivation of Grad-Shafranov and JFKO equations; construction ex- 
act of Euler-potential-based solutions), we rewrite this system of equations in coordinates 
connected with magnetic surfaces. 

Suppose a triply-orthogonal coordinate system (u, v,w) is given, such that the coordinate 
w enumerates magnetic surfaces 2 . Then the pressure P(r) = P(w) is constant on magnetic 
surfaces; the magnetic field has only u— and v— components: 

B = B u e u + B v e v . 

In orthogonal coordinates the metric tensor is diagonal: gij = g%b%y The usual differential 
operators in orthogonal coordinates (u, v, w) have the form 

13/ 19/ 1 df 

gradf =e u ^ = — + e v ^ = — + e w —— — ; (13) 



div A 




^911922933 



(14) 



1 Magnetic surfaces may not exist only in certain cases of field- aligned dynamic MHD equilibria, and in 
the Beltrami case curlB = aB, a = const. A case-by-case classification is found in [H]. 
2 The conditions of existence of such systems are discussed in Remark 1 below. 
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curl A 
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V / fi , Hfi , 225 , 33 



d/du d/dv 



d/dw 

Jg~33~A 3 



(15) 



for any differentiable function f(u, v, w) and vector field A(u, v, w)=Aie u + A 2 e v 
Using the solenoidality condition divB = 0, we get 



d d 

922\/g33B u + TrV9nV933~ B v 



^911922933 \au dv 
which allows to introduce a potential F = F(u, v, w): 
, — / — p dF dF 

Vfi' 22 V5 , 33-D« = — , V5 , llV5 , 33-Dj; = ■ 
OV OU 

Hence the magnetic field is 

1 dF 1 dF 

B = — — e„ — — e v . 

y/922933 dv y/9ll933 °U 

Further, the to-projection of curlB is 0, as both curlB and B are tangent to magnetic 
surfaces w = const. Hence from 1)150: 



0. 



(dF y/gn \ d_ (dF y/gE \ 

V du V '911933/ dv V dv J 9-22933/ 



du V. du y/gTiglxJ dv \ dv ^922933/ 

This is again a potential form: we introduce a potential $ = $(u, v, w) by 

dF J9Y2 9$ dF v^n 

du y/gngw dv ' dv ^922933 du ' 

d 2 F d 2 F 

The compatibility conditions for the derivatives of F must be satisfied: = . 

dudv dvdu 

Hence 

g / g^VgggggP v | d ( d<$> ^gng33 \ = Q 
<9« \<9« v^ii / dv \dv y/g^ J 

The u— and -y— components of the equation curlB x B = /xgradP vanish identically; the 
w— component is 

1 d$ <9 2 $ 1 9$ <9 2 $ 

-/iP (w). (17) 



<7n 9m dudw g 2 2 dv dvdw 

Hence the system of four static isotropic plasma equilibrium equations 1)11)1 ■ in coordinates 
(tt, connected with magnetic surfaces, rewrites as two equations (fT^j). JT7J). 
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Using the definition of gradient and laplacian in non-cartesian coordinates (|13|). the new 
representation can be written as 



A (tt|W )$ = 0, 



<9$ 

?rad M $ ■ grad (Ujl)) — = -fiP (w), (19) 



where (u, v) means that only u— and v— parts of operators are used. 
The magnetic field is expressed through the potential $ as follows: 



1 <9$ 10$ 



B = —p=— e u + -^=— e v = grad (nil)) $, (20) 



and the electric current: 

t 1 m 1 f 1 1 <9 2 $ \ 

J = -curlB = - ^~e u + ^~ e v ■ 21 

A* A* V V922933 ovdw s/gngzz oudw J 

Remark 1. A triply orthogonal coordinate system (u,v,w) with w = const on magnetic 
surfaces can be constructed not for any static equilibrium solution {B,P} of (JTTJ). For a 
family of smooth surfaces w(x, y, z) = const, two other families u(x, y, z) = const, v(x, y, z) = 
const forming a triply orthogonal system can be constructed if and only if the surfaces 
w(x,y,z) = const form a system of Lame 3 [221 123] i-e., the function w(x,y,z) satisfies a 
particular equation of order 3. 

There exist many examples of families of Lame; they include sets of parallel surfaces; sets 
of surfaces of revolution; Ribaucour surfaces, and other families [231- Several appropriate 
examples are discussed below. 

Remark 2. By derivation, every solution {$(«, v, w), P{w)} of the system (|16|) - (J 17)) in 
some orthogonal coordinates (it, v, w) defines a static plasma equilibrium with magnetic field 
flUt satisfying (ITTj) . 

Remark 3. In coordinate systems where gn = gn(u, v), g<n = (u, v), the second equation 
of the system, (JUJ), has a simple energy-connected interpretation. Indeed, the equation can 
be rewritten as 



1 (g rad ( M ,t>)$) 2 = i= p ^ 



w ■ 



v /7^ dw 2 (u ' v > y/g ^ 

3 It is known (result due to Darboux) that for two orthogonal families of surfaces to admit the third one 
orthogonal to both, the two families must intersect in the lines of curvature. A condition that, for a smooth 
family w(x, y, z) = const, there exists a second family of surfaces orthogonal to the given one and intersecting 
it in the lines of curvature, is found in, e.g., |24| . and consists of a PDE of order 3 that w(x, y, z) must satisfy. 
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which is, by (I2(J|) . equivalent to the relation 



9 T+P]=0. (22) 



5f 33 9w V 2 /i 



For incompressible plasma equilibria, the latter means that the component of the gradient 
of total energy density in the direction normal to the magnetic surfaces vanishes. Therefore 
for any MHD equilibrium configuration in which magnetic surfaces w = const form a family 
of Lame, and where gn = gn(u,v), #22 = 5 , 22(w,f), the total energy can be finite only if the 
plasma domain is bounded in the direction transverse to magnetic surfaces. 

For example, plasma equilibria found as cylindrically-symmetric solutions of the Grad- 
Shafranov equation, with domains unbounded in cylindrical radius r and the polar compo- 
nent of B vanishing, are available in literature. For such solutions, in every layer c\ < z < C2, 
the total energy is infinite. However, the magnetic energy in layers c\ < z < c 2 may be finite. 
The same is true for the solutions obtained in, for example, 



Remark 4. As noted by Lundquist [27], the static MHD equilibrium equations (fTTj) are 
equivalent to the time-independent incompressible Euler equations that describe ideal fluid 
equilibria. Therefore static Euler equations may also be presented in the form (J 171) . 

Remark 5. As will be shown in the sections below, in many cases appropriate orthogonal 
coordinates (u, v, w) required by the above theorem may be introduced globally in the plasma 
domain T>. 



3 Exact solutions of Plasma Equilibrium equations in 
the magnetic field-related coordinates and their use 
for modeling 

3.1 Non-trivial Plasma Equilibria arising from in the magnetic 
field-related coordinate representation 

The system of equations (fT^|) . (|17|) under consideration is also essentially non-linear, and 
depends on the metric of unknown orthogonal coordinates. In this section, we list several 
general cases in which explicit solutions of the representation (JTHJ), (JT7J) of Plasma Equilib- 
rium equations can be found. Appropriate examples are found below. 

We are interested in obtaining exact solutions of the system (fTo^) . (fTTj) in different geome- 
tries, but do not restrict ourselves to solutions that can be immediately used as models of 
physical phenomena. Solutions describing force- free (jl2J) and even "vacuum" magnetic fields 
(curlB = 0, P = const) are of interest for physical modeling, because they can serve as 
starting solutions in infinite symmetries and transformations of Plasma Equilibrium equa- 
tions. 

In particular, it was recently shown [SI E] that ideal incompressible MHD equilibrium equa- 
tions possess a Lie group of intrinsic symmetries. If {V(r), B(r), P(r), p(r)} is a solution of 



7 



©j ©5 where the density p(r) is constant on both magnetic field lines and streamlines, 
then {Vi(r), Bi(r), Pi(r), pi(r)} is also a solution, where 



B 1 = 6(r)B + c(r) JjTp V , Y x = B + V , 

Pl = a 2 (r)p, P 1 = CP+ ° 1 . & 2 (r)-c 2 (r) = C = const. 

Here 6(r), c(r) are functions constant on magnetic field lines and plasma streamlines. 
A similar Lie group of transformations exists for incompressible CGL equilibria EI] • 

Another class of transformations worth mentioning here is an infinite-parameter map from 
MHD to CGL equilibrium solutions 0113 EH]. Let {V(r), B(r), P(r), p(r)} be a solution 
of the system (jSJl-d), (jlOft of incompressible MHD equilibrium equations, where the density 
p(r) is constant on both magnetic field lines and plasma streamlines (i.e. on magnetic 
surfaces ^ = const, if they exist.) Then {Vi(r), Bi(r), p±i(r), P||i(r), pi(r)} is a solution to 
incompressible CGL plasma equilibrium system (JBJ-Q, where 

B 1 (r) = /(r)B(r), V^r) = ^(r)V(r), Pl = C oP (v)p/ g 2 (r), 



P±i(r) = Q)/iP(r) + d + (C - / 2 (r)//i) B 2 (r)/2, (24) 



p,n(r) = C pP(v) + C 1 - (Co - / 2 (r)/p) B 2 (r)/2, 

and /(r), g(r) are arbitrary functions constant on the magnetic field lines and streamlines. 
Co, Ci are arbitrary constants. 

For a given equilibrium solution, in the MHD and CGL frameworks, the above transforma- 
tions and the Bogoyavlenskij symmetries can produce a family of solutions connected with it 
but having different behaviour of physical parameters (pressure, density and electric current, 
magnetic and velocity fields). However, such features as solution topology (set of magnetic 
field lines and streamlines), stability with respect to certain classes of instabilities, bounded- 
ness of energy, are inherited by the transformed solutions from the original one [SI 130 EI]- 
If a plasma domain T> G 1R 3 has a boundary, these transformations preserve no-leak-type 
boundary conditions. 

The above-listed transformations can often turn a particular solution into a form suitable 
for physical modeling. However, there exists a lack of diversity of exact equilibrium configu- 
rations that could serve as "starting points" (this lack is observed even in the class of static 
solutions.) Therefore we study the magnetic surface-connected representation of plasma 
equilibrium equations ([lbjh (|17j) aiming at providing methods of building exact solutions in 
different geometries, which could later be transformed into families suitable for modeling. 



3.2 Formulas for exact plasma equilibria in particular geometries 

For a prescribed set of coordinates, the system (fTTj) is a generally nonlinear system 

of equations on two unknown functions v, w), P(w). In coordinates where metric ten- 
sor components are connected in particular ways, formulas defining corresponding plasma 
equilibria can be explicitly written out. 

First we consider cases when v,w) essentially depends on the magnetic surface coor- 
dinate w. Thus the plasma electric current density (|21|) is nonzero. 



Case (A). If the metric tensor components of orthogonal coordinates (u, v,w) satisfy 



911/922 = a{ufb{vfc{wf > 0, g 33 = T' 1 ( w, X(v) - p( u ) ° 2 ^\ r ) , 



where 



then there exists a solution i>, w) of the system 1)1 6|1 . ()17jl in the form 

= C 1 (w)/i(M) + C , 2 HA(t;), (26) 

with the pressure P = const. 

Indeed, under the above assumptions about the relations of u— and v— components of the 
metric tensor, the plasma equilibrium equations ([TB^l. (fTTj) simplify in coordinates (fi,X,w), 
and the solution in the form f)26|) is readily found, together with the necessary expression for 
^33 (& 

The solution defines a force-free plasma equilibrium (|12|) with the proportionality 
coefficient 

1 (2Ci(w) 

a(r) = «(«;) = . 

H w C2\w) dw 

Example. For spherical coordinates (r, 0, <fi) = (w, u, v) the metric g^i — w 2 , (722 = 
w 2 sin 2 m, (733 = 1 satisfies the above relations (J25|) . Hence we find a force- free configuration 

R Ci(r) C 2 (r) dd(r)/rfr 
B = ^— e e + ^— e^, curlB = fiJ = a(r)B, a (r = — — — . 27 
rsinS rsmb 1 ^2( r ) 

with spherical magnetic surfaces r = const. Since the magnetic surfaces are rotationally 
symmetric, several such solutions can be added to produce a non-symmetric force-free plasma 
equilibrium tangent to spheres. 
Moreover, from the equation (JT7j) it follows that for the spherical case 

1 ( fd^V , 1 fd$\ 2 \ 1 r u2 dP(h) JU , 1 



B =— hH +-^ hr = — 2/ ^ — 7r^ d/i + — fl iK w ). 

ur \ \ aw / sin u \ov J I rj dh w 2 



where a\(u, v) is generally not identically zero, and is never identically zero for force-free 
plasmas (since B 2 does not vanish in the whole plasma domain). Therefore any force-free 
and general non-force- free plasma equilibrium configurations with spherical magnetic surfaces 
have a pole-type singularity at the origin r = 0, and infinite magnetic energy, if the plasma 
region includes the origin. 

Case (B). If the metric tensor components of orthogonal coordinates (u,v,w) satisfy 

gn /g 22 = a(u) 2 b(v) 2 , g 33 = F 2 (w), (28) 

then 

$!(«, v, w)= I t(k) [dHe™^ cos(n(A;)A) + C 2 (w)e Tn{k)l1 sin(n(fc)A)] dk, (29) 



$ 2 (n, v, w) = t{k) [C 1 (w)e rn{k)x cos(n(fc)//) + C 2 (w)e Tn{k)x sin(n(A;)^)] dk, (30) 



define solutions to the system of isotropic plasma equilibrium system (j!6|) . (fTTj) . for any 

CAw) and CAw) satisfying — (C?(w) + C^iw)) = 0. Here again fi(u) = f a(u)du, \(v) = 

dw 

J l/b(u)du, r = ±1, n(k) is an arbitrary function, and t(k) is an arbitrary generalized 
function (for each solution, n(k),t(k) must be chosen so that the integral converges). 
These solutions also correspond to force-free plasma equilibria (|12j) with the coefficient 

1 dCi{w) 

a(r) = a(w) 



J r (w)C2(w) dw 



This statement follows from noticing that the relations (|28|) turn the plasma equilibrium 
system flTB|) . (|T7|) into a usual Laplace equation and a total derivative with respect to Wi = 
f J r {w)dw\ 

&w + $aa = 0, + $ a $a Wi = 0; 

the solutions flU, © follow. 

Remark 1. Relations of the type (J23j) and (|2*H|) between metric components coefficients are 
not unnatural. The simplest example is coordinates obtained by a conformal transformation 
of the complex plane, as shown in Section [Ql below. 

We now turn attention to the remaining case when the unknown function Q(u,v,w) is 
independent on the magnetic surface variable w. 

Case (C). As noted above, the first plasma equilibrium equation (fT7)|) in magnetic surface 
coordinate representation coincides with the (u — v)— part of a Laplacian (|18|). Hence in any 
coordinate system where the 3D Laplace equation A( U)V)W )<p(u,v,w) = admits a solution 
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independent of one of the variables (w), there exists a gradient ("vacuum") magnetic field 
configuration 



divB = 0, curlB = 



(31) 



corresponding to this solution, and this magnetic field is tangent to surfaces w = const. 

Remark 2. Availability of "vacuum" magnetic fields. Many classical and esoteric 
coordinate systems admit geometrically nontrivial two-dimensional solutions of the Laplace 
equation, as found in literature, for example, (2H1- New systems of coordinates may be 
constructed where the Laplace's equation will be separable or have two-dimensional solutions. 
The list of conditions on the metric coefficients necessary and sufficient for separability of 
Laplace equation and existence of two-dimensional solutions is available in j2H]. In the same 
book one finds methods of producing new triply orthogonal coordinate systems by conformal 
transformations of the complex plane. 

Remark 3. Use of "vacuum" magnetic fields. Magnetic fields of the type (|31j) can 
be found in different geometries. Under the action of Bogoyavlenskij symmetries (|23jl or 
transformations to anisotropic CGL equilibria (|24J1. these gradient fields are transformed into 
ones with curlB ^ and give rise to non-trivial dynamic MHD and static and dynamic CGL 
plasma equilibria in the same geometry as the original "vacuum" solution. The examples 
are given below. 

3.3 Construction of additional plasma equilibria using the mag- 
netic surface-related coordinate representation 

The following two statements extend classes of solutions of static MHD equilibrium equations 
(|16 [) -([17 )) in magnetic-surface-related coordinate representation. Namely, under particular 
conditions on the metric, the equilibrium magnetic field can be extended with a Killing 
component in the w— direction. 

The first statement presents a transformation of a "vacuum" curl-free magnetic vector field 
depending only on two variables into an extended magnetic vector field which is not force- 
free or gradient, and thus gives rise to a new non- degenerate solution to plasma equilibrium 
equations The example of use of this transformation for the extention of a solution 

class is found in the subsequent section |4~B1 

Statement 1 If <p(u,v) is a solution to the system Mb)) ,\ll\j in coordinates (u,v,w) with 



then not only the magnetic field A2U\) with pressure P(w) = const solves the Plasma Equilib- 
rium equations 177)) . but so does the extended magnetic field 



properties 



gil = gil(u,v) = # 22 , 033=033(w), 



(32) 



B 




1 d<f) 



1 dcj) 



(33) 



w 
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with pressure 



p = c - K\u,v)/2, (34) 

where K(u, v) satisfies 

d 2 K{u,v) d 2 K{u,v) n JJL , , . . 

— ^2 1 ^2 = °' S rad 0( M > w ) ' & rad u ) = °- 

This statement is verified directly by substituting the magnetic field (JHSJ) and the pressure 
(I34j) into the static plasma equilibrium system The function K(u,v) is a harmonic 

conjugate of the solution (f>(u,v). 

It turns out to be possible to add a w— component to a wider class of solutions of plasma 
equilibrium equations ()16 |) -([T7 |) in magnetic-surface-coordinate representation. The following 
statement extends a "vacuum" curl-free magnetic vector field depending only on two variables 
and tangent to surfaces w = const to "vacuum" fields that have non-zero w-components. 
The conditions on the metric in this case are more relaxed. 

Statement 2 If<f)(u, v) is a 2- dimensional solution to the plasma equilibrium system ilb}) , \lT\ i 
in the coordinates (u,v,w) with properties 

gu=gu(u,v), #22 = gw(u,v), g 33 = a(w) 2 J r2 (n, w), (35) 

then not only the magnetic field \2(ty with pressure P(w) = const solves the Plasma Equilib- 
rium equations ill]) , but so does the extended magnetic field 

1 dcj) 1 dcj) D 
B= H-^ + ir u ^ + H- w e ^ D = ™™t. (36) 

This magnetic filed is a vacuum magnetic field: divB = 0, curlB = and corresponds to 
plasma equilibria with P = const. 

This transformation is used in the example in Section l4~2*l below. 



4 Examples of exact plasma equilibria. Physical mod- 
els. 

4.1 Isotropic and anisotropic plasma equilibria tangent to nested 
ellipsoids 

In this example, we construct a family of generally non-symmetric plasma equilibria with 
ellipsoidal magnetic surfaces. 

We start from the construction of vacuum magnetic fields tangent to ellipsoids, using the 
magnetic-surface-connected representation of plasma equilibria equations (fTB^l . (JT7J) (Section 
13.21 Case C.) Then transformations are applied to this gradient solution to produce families 
of non-trivial isotropic and anisotropic plasma equilibria. 
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An application of the resulting solutions to modeling solar photosphere plasma near active 
regions is discussed. 

(i). A vacuum magnetic field tangent to ellipsoids. The ellipsoidal coordinates are 



u = e, b 2 <9 2 < c 2 , 

v = A, < A 2 < b 2 , 

W — T), C 2 < T] 2 < + OO. 

The coordinate surfaces are 

= 1, (ellipsoids,?] = const), 

= 1, (one-sheet hyperboloids, 9 = const), 

= 1, (two-sheet hyperboloids, A = const). 

Laplace's equation is separable in ellipsoidal coordinates, and we take a solution depending 
only on (9, A), so that its gradient has zero ^-projection transverse to ellipsoids, but is tangent 
to them: 



X 2 


+ 


y 2 




z 2 


V 2 
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-c 2 
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Here sn(x, k) is the Jacobi elliptic sine function. The inverse of it is an incomplete elliptic 
integral 

f z 1 

Feii(z,k)= / = dt. 

1 ; Jo VT^Vi - k 2 t 2 

«3>i(#, A) does not depend on w, and therefore evidently satisfies both equations (fTHjl . (fT7|l . 
The resulting magnetic field (j^Oj) is tangent to ellipsoids 77 = const, and has a singularity at 
9 = A, i.e. on the plane y — 0. 

However one may verify that for a plasma region c < rji < t] < ?] 2 the total magnetic energy 
f v B 2 /2dv is finite. Also, if one restricts to a half-space y > or y < 0, then the magnetic 
field is well-defined in a continuous and differentiable way. 

If the magnetic field is tangent to the boundary of a domain, one can safely assume that 
outside of it B = identically. This is achieved, as usual, by introducing a boundary surface 
current 

i 6 (ri) = fi -1 B(ri) x rw(rx), (37) 

where ri is a point on the boundary of the domain, and n out is an outward normal. 

Fig. Q shows several magnetic field lines for the case (b = 7, c = 10, A\ = A2 — 0, Bi = 
1/100, B 2 = 1/30) on the ellipsoid rj = 12. For this set of constants, the vector of the 



13 



magnetic field has the form is 




'100-02 

Bo = , =e e - - V - =^e A (3* 

y/(0 2 - \ 2 )( V 2 -8 2 ) s/(v 2 ~ A 2 )(0 2 -A 2 ) 

This "vacuum" (gradient) magnetic field is used to produce non-trivial dynamic isotropic 
and anisotropic plasma equilibria, as shown below. 

We remark that though the magnetic field lines of the field ()38|) have a plane of symmetry 
x = 0, a non-zero choice of constants (Ai, A 2 , B>i, P> 2 ) would produce a completely non- 
symmetric magnetic field tangent to a family of ellipsoids. 

(ii). Isotropic dynamic plasma equilibrium with ellipsoidal magnetic surfaces. 

The above vacuum magnetic field B is indeed a trivial solution to the general isotropic 
plasma equilibrium system ©-© with V = 0, P = Pq = const and an arbitrary density 
function p = po(r). 

If we choose po( r ) to be constant on magnetic field lines (plasma streamlines do not exist as 
there is no flow), then the infinite-parameter transformations (|23|) become applicable to such 
configuration. Applying them formally, we obtain a family of isotropic plasma equilibria 

n( r l 

Bi = m(r)B 0) V x = — -±L==B , 

a(r)y/fipo(r) 



Pl = a 2 (r)p (r), P 1 = CP - n 2 (r)B 2 /(2fj,), (39) 



m 2 (r) — n 2 {r) = C = const, 

where a(r), m(r), n(r), Po( r ) are functions constant on magnetic field lines and streamlines 
(which coincide in this CcLSC, clS Vi and B x are collinear). 

We consider a plasma configuration in a region T> in the half-space y > between two 
ellipsoid shells rjx, rj 2 : c < r)\ < rj < r/ 2 , (we take c = 10; c is one of the parameters of the 
elliptic coordinate systems used for this solution.) 

Outside of the region, we assume Bi = 0, by introducing a corresponding surface current 
(|37|). We also assume Vi = 0, which can be done because the streamlines are tangent to the 
boundary of the plasma domain T>. 

The magnetic field lines in the chosen region are not dense on any 2D surface or in any 3D 
domain, therefore the arbitrary functions a(r), m(r), n(r), Po( r ) can be chosen (in a smooth 
way) to have a constant value on each magnetic field line, thus being in fact functions of two 
variables enumerating all magnetic field lines in the region of interest (for example, rj and A, 
which specify the beginning of every magnetic field line). 
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Figure 1: A magnetic field tangent to ellipsoids. 

Lines of the magnetic field (|38[) tangent to the ellipsoid t] = 12. The shown ellipsoid is 
a magnetic surface from the family of nested ellipsoids rj = const in classical ellipsoidal 
coordinates. 

This configuration is smoothly defined in the half-space y > 0. 
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We remark that unlike the initial field B , the vector fields Bi and Vi are neither potential 
nor force-free: for example, curlBx = gradm(r) x B Bi. But both Bi and Vi satisfy the 
solenoidality requirement. 



Direct verification shows that, with a non-singular choice of the arbitrary functions, the 
total magnetic energy E m = 1/2 J y B\ dv and the kinetic energy E k = 1/2 f v p V 2 dv are 
finite. The magnetic field, velocity, pressure and density (Bi, Vi, pi, Pi) are defined in a 
continuous and differentiable way. 

The presented model is not unstable according to the known sufficient instability condition 
for incompressible plasma equilibria with flows proven in [3U] (see also jH].) The latter states 
that if V B, then a plasma equilibrium with constat density is unstable. In the presented 
example, V || B (the density pi can be chosen constant). 

(iii). Anisotropic plasma equilibrium with ellipsoidal magnetic surfaces. When 
the mean free path for particle collisions is long compared to Larmor radius, (e.g. in strongly 
magnetized plasmas), the tensor-pressure CGL approximation should be used. The model 
suggested here describes a rarefied plasma behaviour in a strong magnetic field looping out 
of the star surface. 

To construct an anisotropic CGL extension of the above isotropic model, we use the 
transformations (|24|) (Chapter 3) from MHD to CGL equilibrium configurations. Given 
Bi, Vi,Pi,pi determined by (J39|) with some choice of the arbitrary functions a(r), m(r), 
n(r), po( r ); we obtain an anisotropic equilibrium B 2 , V 2 ,p|| 2 , p_|_ 2 , P2 defined as 



/(r), g(r) are arbitrary functions constant on the magnetic field lines and streamlines, i.e. 
again on constant on every plasma magnetic field line, and Co, C\ are arbitrary constants. 
Setting P = in (|3Uj) and making an explicit substitution, we get 



B 2 = /(r)B 1; V 2 = y(r)Vi, p 2 = C oP ip/g{r)\ 



VL2 = C pPi +d+ (C - /(r) 2 //i) Bj/2, 



(40) 



pp = CopPi + Ci- (Co - f(r) 2 /p) Bj/2 



B 2 



/(r)m(r)B , V 2 = g(r) 



n(r) 



B , P2 = C a 2 (v)p (r)p/g(r) 2 , 



a 



(r)v^(r) 



RL2 



C 1 + ^(C Cp-/ 2 (r)m 2 (r)) 



(41) 



2 



Ci + -° (/ 2 (r)m 2 (r) - C C/i - 2C n 2 (r)) . 
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It is known ([TSJ E3) that for the new equilibrium to be free from a fire-hose instability, 
the transformations (|24j) must have Co > 0. 

is the pressure component perpendicular to magnetic field lines. It is due to the rotation 
of particles in the magnetic field. Therefore in strongly magnetized or rarified plasmas, where 
the CGL equilibrium model is applicable, the behaviour of p± should reflect that of B 2 . 

In the studies of the solar wind flow in the Earth magnetosheath, the relation 

P±/P\\ = 1 + 0.847(5 2 /(2p||)) (42) 

was proposed [HI]. We denote k(r) = CoC/i — / 2 (r)m 2 (r) and select the constants and 
functions Co, C, /(r), m(r) so that k(r) > in the space region under consideration. From 
(jUJ, we have: 

P± 2 -P|| 2 = |j(2Mr) + 2Con 2 (r)), 

or 

P±2 =1 + 2fc(r)+2C n 2 (r) B 2 
Pp /i/ 2 (r)m 2 (r) 2pp' 

which generalizes and includes the experimental result (J42I) . 

(iv). A model of plasma behaviour in arcade solar flares. 

Solar flares are known to take place in the photospheric region of the solar atmosphere and 
are connected with a sudden release of huge energies (typically 10 22 — 10 25 J) (e.g. |H2j, PP- 
331-348). Particle velocities connected with this phenomenon (about 10 3 m/s) are rather 
small compared to typical coronal velocities (~ 5 ■ 10 5 m/s), therefore equilibrium models 
are applicable. 

Morphologically two types of solar flares are distinguished: loop arcades (magnetic flux 
tubes) and two-ribbon flares. Flares themselves and post-flare loops are grounded in from 
active photospheric regions. 

As noted in [32], p. 332, "rigorous theoretical modelling has mainly been restricted to 
symmetric configurations, cylindrical models of coronal loops and two-dimensional arcades. " 

The configurations described in (ii) and (iii) can serve as non-symmetric 3D isotropic and 
anisotropic models of quasi-equilibrium plasma in flare and post-flare loops, where magnetic 
field and inertia terms prevail upon the gravitation potential term in the plasma equilibrium 
equations: 

V x curl V ^> grad<p, — B x curl B ^> pgrady?. 

H 

where </? is the star gravitation field potential. 

The relative position and form of the magnetic field lines in the model, with respect to 
the star surface, are shown on Fig. El The characteristic shape of the magnetic field energy 
density B and the pressure P along a particular magnetic field line, for the isotropic case 
(ii), are given on Fig. El 
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Figure 2: A solar flare model - magnetic field lines. 



The model of a solar flare as a coronal plasma loop near an active photospheric region. The 
position and shape of several magnetic field lines are shown with respect to the star surface; 
magnetic field is tangent to nested half-ellipsoids. 
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Figure 3: A solar flare model - plasma parameter profiles. 

A model of a solar flare - a coronal plasma loop near an active photospheric region. 

The figure shows the characteristic shape of the magnetic field energy density B and the 

pressure P curves along a particular magnetic field line. (The isotropic case.) 
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Magnetic field lines in the model are not closed; therefore by introducing a surface current 
of the type (J37J), a plasma domain V can indeed be restricted to any flux tube, with boundary 
tangent to magnetic field lines, and the magnetic field can be chosen zero outside (together 
with the velocity in models (ii), (iii)) by the introduction of a boundary surface current. 
The current sheet introduction is not artificial - as argued in [33], in a general 3D coronal 
configurations the current sheets between flux tubes are formed (see also: [221, p. 343.) 

The isotropic MHD model (ii) is valid when the mean free path of plasma particles is much 
less than the typical scale of the problem, so that the picture is maintained nearly isotropic 
via frequent collisions. 

However, the CGL framework must be adopted when plasma is rarefied or strongly mag- 
netized. For such plasmas, we propose the anisotropic model (iii), for which the requirement 
of plasma being rarefied can be satisfied by choosing a(r) sufficiently small. 

4.2 Isotropic and anisotropic plasma equilibria in prolate spheroidal 
coordinates. A model of mass exchange between two distant 
spheroidal objects 

In this example, families of non-symmetric exact plasma equilibria is prolate spheroidal co- 
ordinates with finite magnetic energy are obtained, in isotropic and anisotropic frameworks. 
On the basis of these solutions, a model of the quasi-equilibrium stage of mass exchange by 
a plasma jet between two distant spheroidal objects is suggested. 

(i). Vacuum magnetic configuration in prolate spheroidal coordinates. Consider 
the prolate spheroidal system of orthogonal coordinates: 

U = 9, < 9 < 7T, 

v — (f), < (f) < 2ir, 
w — T), < rj < +oo. 

The coordinate surfaces are 

(prolate spheroids,?] = const), 

(two-sheet hyperboloids, 9 = const), 
(half planes, = const), 

and the metric coefficients 

9rp) = gee = a 2 (sinh 2 rj + sin 2 9), g^ = a 2 sinh 2 1] sin 2 9. 

It is known that the 3-dimensional Laplace equation is separable in this system, and it 
admits the axially-symmetric family of solutions of this equation [2*%] : 

$ = H(r])T(9), 



x 



a 2 sinh 2 rj a 2 sinh 2 77 a 2 cosh 2 77 



x 



y 



+ 



a 2 sin 2 9y a 2 sin 2 9 a 2 cos 2 9 

tan(0) = — , 

x 
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H(rj) = AiP p (coshr/) + B^coshr]), T(0) = A 2 V p (cos9) + B 2 C p (cos9), 

where C p (z) = C p (z), jC p (z) = C p (z) are theLegendre wave functions of first and second 
kind respectively. 

This solution evidently satisfies both plasma equilibrium equations in magnetic-surface- 
related coordinates (|Tfi|) . (fTTjl : the first one because the usual and truncated Laplace equa- 
tions coincide when $ is a function of two variables, and the second - identically due to the 
independence of $ on w. 

In the case of integer p, the above 2-dimensional solution expresses in ordinary Legendre 
functions of the first and second kind: 

H( V ) = AxPpicoshr]) + B^coshr]), T{9) = A 2 P p (cos6) + B 2 Q p (cos6). (43) 

From the above family, a particular axially symmetric function $(77, 9) with an asymptotic 
condition 

lim $(77, 0) = M z, (44) 

|r|— >oo 

can be chosen; its gradient is asymptotically a constant vector field in the cartesian z- 
direction: grad$(77, 6) = M e z . 
This solution has the form 



1 (cosh 77) I 

n cosh 77 — cosh?7o— — : r > , 

Qi (cosh 770) J 

and the corresponding magnetic field B = grad<3>(7?, 8) is 

B = grad$ (v,0) = — ° W ' J e n + — °^ } eg. 45 



The magnetic surfaces this field is tangent to are nested widening circular tubes along z- 
axis, perpendicular to the spheroid rj = r] and asymptotically approaching circular cylinders 
x 2 + y 2 = const. A graph for the choice {a = 2, M = 1, 770 = 0.3} with a single magnetic 
field line shown is presented on Fig. 0] 

Each tube is uniquely defined by the value of 9 of its intersection with the base spheroid. 
The two shown on Fig. 0] correspond to 61 = 0.07 and 9 2 = 0.12. 

If a "winding" polar component 
D 



a sinh 77 sin 6, D = const, 



is added to the field (J43j) . which can be done by Statement El of section 13.31 above, a new 
vacuum magnetic field is obtained: 

B w = grad $ 77, 0) = = y + — 01 ' e e + -— e*. (46) 
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Figure 4: A magnetic field flux tube normal to a prolate spheroid. 

The magnetic field (|45jl and magnetic surfaces in prolate spheroidal coordinates. The mag- 
netic surfaces are nested widening circular tubes along the z-axis, perpendicular to the 
spheroid 77 = 770 and asymptotically approaching circular cylinders x 2 + y 2 = const. 
The graph is built for the choice {a = 2, Mq = 1, 770 = 0.3}. The two magnetic surfaces 
shown here correspond to 6\ = 0.07 and 6 2 = 0.12. 
A sample magnetic field line on the outer surface (8 2 = 0.12) is plotted. 
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Figure 5: A winding magnetic field in prolate spheroidal coordinates. 

Two opposite field lines of the magnetic field (|4fi|) winding around a magnetic surface are 
shown; a solution is constructed in prolate spheroidal coordinates. 

The magnetic surfaces are nested widening circular tubes along the z-axis, perpendicular to 
the spheroid rj — rjo and asymptotically approaching circular cylinders x 2 + y 2 = const. The 
magnetic surface shown here is defined by its intersection with the spheroid at 9 = 0.3. 
The graph is built for the choice {a = 2, M = 1, r] = 0.3}, D = 2.13. 



23 



A graph for D = 2.13 showing two opposite field lines winding around a magnetic surface 
is presented on Fig. 

The physical model. We use the vacuum solutions with and without the polar component, 
(|45[) and (J46|) . to model a quasi-equilibrium process of mass exchange by a plasma jet between 
two distant spheroidal objects. 

The useful property of the solutions is that their magnetic surfaces tend to cylinders by 
construction. Also, if B is a vacuum magnetic field, then (— B) is a vacuum magnetic field, 
too. 

Hence one may effectively glue one copy of such solution with another copy, the latter 
being rotated on the angle tt with respect to an axis orthogonal to the axis of the magnetic 
surface, translated on the distance much longer than the size of the initial spheroid (Fig. |)J), 
and taken with the opposite sign. 

It is possible to show that for a given solution ( (J45|) or (|46|) ) . the rate of growth of the 
tube radius B r /B z has the leading term z~ 3 at z — > oo, hence the magnetic field lines of the 
"glued" solution will not have significant "cusps" - discontinuities of derivatives. 

The resulting force-free vacuum magnetic field can be used to construct isotropic plasma 
equilibria with flow by virtue of the Bogoyavlenskij symmetries (|23|). or anisotropic plasma 
equilibria with and without flow, with the help of the MHD— >CGL transformations (|24j). 

For example, after the application of the Bogoyavlenskij symmetries to the field ~B W , one 
gets an isotropic dynamic configuration 

Bx = m(r)B„ Vi = —-^d==B w , 

Pl = a 2 (r)p (r), P 1 = CP - n 2 (r)B 2 / (2//). (47) 



m 2 (r) — n 2 {r) = C = const, 

where a(r), m(r), n(r), Po( r ) are functions constant on magnetic field lines and streamlines 
(which coincide, as Vi and Bx are collinear). 

To construct an anisotropic CGL extension of the above isotropic model, we again use 
the transformations (J23|. The resulting anisotropic equilibrium B 2 , V 2 ,P||2,Ri_2, P2 is then 
defined by (P was set to 0): 

B 2 = /(r)m(r)B tt , V 2 = g(r)—-±L=B w , p 2 = C a 2 (r)p (r)^/g(r) 2 , 

«( r )V/ i Po(r) 

p±2 = C x + ^ (CoC/i - f(r)m 2 (r)) , (48) 
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Figure 6: A model of mass exchange between two distant spheroidal objects by a plasma jet. 

The magnetic surfaces are nested widening circular tubes along the z-axis, perpendicular 
to the spheroid f] = f]o and asymptotically approaching circular cylinders x 2 + y 2 = const: 
B r /B z = 0(z~ 3 ) at z -* oo. 

Shown here is the procedure of gluing one copy of a solution ( (J46|) or ()45Jl) with another 
copy, rotated on the angle ir with respect to an axis orthogonal to the axis of the magnetic 
surface, translated on the distance much longer than the size of the initial spheroid, and 
taken with the opposite sign. 
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n2 = CJl + lt ( /2(r)m2(r) " c ° Gfi ~ 2c ° n2 ^) ■ 

The physical requirements and applicability bounds are the same as described in the pre- 
vious model (see sec. 14.1)1 . The relation between the pressure components of anisotropic 
pressure tensor, is also the same: 

P±2 2fc(r) + 2C n 2 (r) B 2 , 2 2 

= 1 + Wi \ 2f \ o — ' fc ( r ) = c oC^-f (r)m (r), 

P|| 2 /i/ 2 (r)m 2 (r) 2p|| 2 

which is in the agreement with the observation-based empiric formula (|42jl. 

We remark that the same way as in the previous model, the values of all the arbitrary 
functions of the transformations f)23jl. ()24)1 can be chosen separately not on every magnetic 
surface, but on every magnetic field line. Thus these free functions are actually functions 
of two independent variables specifying the origin of every magnetic line on the starting 
spheroid, and the resulting exact solution has no geometrical symmetries. 

If the constant D in the initial field B w (J46)) is different from zero, then the w - component 
of this field has a singularity on the z-axis, and the plasma domain T> must be restricted 
to a volume between two nested magnetic surfaces so that the z-axis is excluded (see Fig. |7J). 
However, the families of transformed isotropic (|4"7j) and anisotropic (J48)) magnetic fields Bi, 
B2 are smooth everywhere, if the non-singular field Bo (145)1 is used instead of B„. Then the 
plasma domain T> can be chosen to be a region inside any flux tube or between two nested 
ones (Fig. 0) The domain on Fig. EH is simply connected; the one on Fig. 0d is not simply 
connected. The axis of symmetry of the magnetic surfaces coincides with the big axis of 
both spheroids. 

In dynamic isotropic and anisotropic cases, by the properties of solutions constructed from 
static configurations V = by the transformations (|23)) or 1)24)1 . the plasma velocity has 
the same direction as the magnetic field, V || B, so the configuration is interpreted as a 
magnetically driven matter flow from one spheroid to another. 

The presented solution models the quasi-equilibrium stage during the time interval T, with 
the requirement 

T ■ S max ■ max \pV\ < M , 

where M is the mass of the spheroid objects, and S max is the area of the maximal section 
of the plasma domain transverse to the flow lines. 

4.3 Example 3. Generation of orthogonal coordinate systems by 
coordinate transformations 

Given a coordinate system (x 1 ,^ 2 ,^ 3 ) in the flat space 1R 3 that satisfies Riemann equations 
Rijki — 0, one can use an arbitrary coordinate transformation 

u l = u i (x\x 2 ,x 3 ), 2 = 1,2,3, 
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Figure 7: A model of mass exchange between two distant spheroidal objects by a plasma jet: 
possible plasma domains. 

Examples of possible plasma domains T> for isotropic and anisotropic equilibria models (sec. 

El) 
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and, by tensor transformation rules, the Riemann tensor of the resulting coordinates will 
also be identically zero: 

Qoc^ 1 c)oc^ Q^jC^ 
R'ijki = Rabcd -^Q^Q^kQ^I = °> i,3,k,l,a,b,c,d= 1,2,3. 

Here we give an example of transformations that produce orthogonal coordinates and have 
metric coefficients satisfying the sufficient condition for a force-free plasma equilibrium of 
the type (J2BJ to exist. 

Consider the plane transformations 

x = £ 1 (u,v), j/ = 6j(u,v) (49) 
satisfying the Cauchy-Riemann conditions 

dx dy dx dy , . 

du dv' dv du 

(Here x, y are cartesian and u, v curvilinear coordinates). The property of conformal map- 
pings is that it preserves angles, hence the families of curves u = const, v = const in the 
plane are mutually orthogonal. 
If we consider the corresponding 3D cylindrical mapping 

x = £i(u,v), y = &(u,v), z = w, (51) 

it defines an orthogonal coordinate system with metric coefficients 

911 = 922 = q + k ' ^33 = 1- (52) 



du J V du 



First type of solutions. These metric coefficients exactly satisfy the conditions used in 
Cases (A) and (B) of Section I3~2*l hence and in the coordinates (u, v, w) a force- free magnetic 
field exists: 

Many examples of such cylindrical transformations can be suggested. The simplest ones 
include power, logarithmic, exponential, hyperbolic, elliptic and other types of conformal 
complex plane mappings. We do not consider solutions of this type in detail here. 

Remark. The field lines of force-free magnetic fields (|53jl lie in planes z = const. A constant 
z- component can be added to the fields of this type as follows: 

B= fC^C^ \ d { ct( W ) + Cl(w)) = 0, D = const. (54) 

Then the electric current density J = curlB//^ does not change, and the equilibrium remains 
force- free. 
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Second type of solutions. The solutions presented below can be built in any coordinate 
system obtained from the cartesian coordinates (x, y, z) by a conformal plane transformation 

In the transformed coordinates (|51j). the metric coefficients are (|52j). hence the complete 
Laplace equation in coordinates (u,v,w) evidently has 2-dimensional solutions <$>(u,v): 

<9 2 $ <9 2 $ 

du 2 dv 2 

Therefore, by Case (C) of Section l3~2*l in these coordinates there exists a vacuum magnetic 
field (J20| in the z = const-plane defined by 
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which can be given a non-trivial z-component using Statement Q i n Section 13.31 

/ 1 <9$ 1 (9$ \ 

B = — — — . — — . K(u,v)\, AK(u,v)—0, grad^fw, v) -graAKiu, v) — 0. 

V V9n ou Jg^ dv J 



By Lemma El such field and the pressure P(u,v) = C — K 2 (u,v)/2 (C = const) satisfy 
the full plasma equilibrium system (fTTj) 

curlB x B = ugradP, divB = 0. 



Example. We now give a particular example in elliptic cylindrical coordinate system defined 
by a conformal transformation that acts on the complex plane as Z' = a cosh Z: 

x = a cosh u cos v, y = a sinh u sin v, z — w. 
We let a = 1 and choose a function satisfying A( Ujl ,)$(u, v) = 0: 

u) = sinh u cos u + 0.1 sinh 2u cos2w — 3v + Ci. 
A conjugate harmonic function for it is 

K(u,v) = cosh u sin v + 0.1 cosh 2u sin 2-u + 3u. 

The level curves K(u,v) = const are presented on Fig. |H]and coincide with the projections 
of the magnetic field lines on the (x, |/)-plane. 

The corresponding plasma equilibrium solution on the cylinders K(u, v) = const has a 
simple representation 

_ cosh u cos v + 0.2 cosh 1u cos 2v — sinh u sinv — 0.2 sinh 2m sin 2v — 3 

= e " ^ 1 Gv 

V cosh 2 u — cos 2 f v cosh 2 it — cos 2 i> (55) 

+ (3w + cosh m sin t> + 0.1 cosh 2m sin 2t>)e z 
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Figure 8: An example of non-circular cylindrical magnetic surfaces. 



A family of cylinders with non-symmetric closed cross-sections that are the magnetic 
of sample isotropic and anisotropic plasma equilibria found in sec. 14.31 
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The u- and f-components of this magnetic field evidently have a singularity at u = v = of 
the order p" 1 , where p = \/u 2 + v 2 is the "distance" to singularity. In cartesian coordinates, 
the singularity is located at (x — ±1, y — 0). 

The levels of the magnetic surface function K(u, v) = const encircle the singularities, and 
degenerate into a line segment as K(u, v) — > 0, as seen from Fig. |H] 

If the magnetic field with pressure P(u, v) = C — K 2 (u, v)/2 (J3"4"|) is to be used in any 
isotropic model, the plasma domain T> is to be restricted to the cylindrical volume between 
any two magnetic surfaces: V = {(u,v) : < K\ < K(u,v) < K 2 . Using the fact that the 
magnetic field is tangent to the surfaces K(u,v) = const, this can be done by introducing a 
boundary surface current (|37|). Then outside of T> magnetic field is zero: B = 0. 

However, if one is to use the described configuration as an initial solution for building an 
anisotropic static equilibrium by virtue of the transformations (|24|) . then he should select 
the function /(r) (that must be constant on magnetic surfaces, /(r) = f(K)) as follows: 

f(K) = ±——a(K), 

max { \B 1 , \B 2 } 

K{u,v)=K 

where a(K) is some function with compact support and the property \ k=o = 0, and 
\Bi\, \B 2 \ are respectively the u- and w-components of the field ()55|). 

This function is evidently finite in any domain T> bounded by a level K(u, v) < K. On the 
degenerate surface K(u, v) — by continuity we have B = 0. On the outer boundary of the 
domain, a surface current (|3*Tj) must be introduced, to ensure so that B = outside of T>. 

Remark. Vacuum magnetic fields in rotational coordinate systems. 

Among the classical and esoteric coordinate systems where the Laplace equation is sepa- 
rable or R-separable [2B1> many are rotationally symmetric systems, with metric coefficients 
independent of the polar angle <fi. 

In all such systems, the Laplace equation has solutions independent of <fi. Examples are 
toroidal coordinates, usual and inverse prolate and oblate spheroidal coordinates, cap-cyclide, 
disk-cyclide, cardioid coordinates and several others (see |28j). 

By the statement of Case (C) of Section l3~2l gradient "vacuum" magnetic fields can be built 
in such coordinates, tangent to the magnetic surfaces, which are in this case vertical half- 
planes = const. Again, a non-zero 0-component can be added to these vacuum magnetic 
fields, to make them non-planar. 

Examples of such vacuum magnetic fields and the corresponding isotropic and anisotropic 
plasma configurations obtained from them by transformations ()23|) . ()24|) will be built in 
consequent papers. 

Particular solutions obtained in different coordinates can have simple algebraic representa- 
tion only in the corresponding coordinates, therefore different rotational coordinate systems 
may not be considered equivalent, from the computational point of view. 

We also remark that the magnetic fields constructed this way can not be found from Grad- 
Shafranov equation, which describes plasma equilibria with magnetic surfaces \l/(r, z) = 
const, whereas in the above case magnetic surfaces are ^ = ^(<p) = const. 
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5 Conclusion 



In this paper, a method of construction of exact plasma equilibria is presented. It is used 
for producing dynamic and static equilibria in different geometries, in both classical MHD 
and anisotropic tensor-pressure CGL frameworks. 

The method is based on representing the system of static isotropic plasma equilibrium 
equations (jllj) in coordinates (u,v,w), such that magnetic surfaces coincide with the coor- 
dinate level surfaces w = const. Such representation is valid when the family of magnetic 
surfaces can be a part of a triply orthogonal system of surfaces, i.e. when it forms a family 
of Lame. In such coordinates, the system of four static Plasma Equilibrium equations (fTTf 
is reduced to two partial differential equations for two unknown functions. The first of the 
equations of the system is a "truncated" Laplace equation (J16)) . and the second one, (J17|) . 
has an energy-connected interpretation (Section |21) 

Instead of four unknown functions of the static MHD equilibrium system, B(r) and -P(r), 
that depend on three spatial variables, the new system of equations employs only two func- 
tions - v,w) and P(w), and the magnetic field B(r) is reconstructed from the relation 
flUJ). 

In Section |3J sufficient conditions on the metric coefficients are established under which 
exact solutions of particular types can be found in corresponding coordinates. In particular, 
if the conditions (|23j) or (|2*Hj) are satisfied, a force- free plasma equilibrium {P{w) = const) 
can be constructed. 

In coordinates where the general 3D Laplace equation admits 2-dimensional solutions, 
"vacuum" magnetic fields divB = 0, curlB = tangent to magnetic surfaces w = const) 
can be built; such field often have non-trivial geometry. Statements^ and El extend this class 
of solutions in some coordinate systems. 

The use of "vacuum" gradient fields for plasma equilibrium modeling is discussed. Such 
solutions can serve as initial solutions in the infinite-parameter Bogoyavlenskij symmetries 
(|2*3*|) . or in the MHD — > CGL transformations (j21J). The application of these symmetries 
and transformations, for each initial solution, generates families of dynamic MHD and static 
and dynamic CGL equilibria; suitable solutions from these families can be used as physical 
models. 

We use the described procedure to construct particular examples of plasma equilibria. The 
first example (Section 13 .2|) is a set of non-Beltrami Force-Free plasma equilibria ()12|) in a 
prescribed geometry - with spherical magnetic surfaces and the force-free coefficient a(r) 
(TH?1) being a function of the spherical radius. When used in a linear combination, these 
solutions give rise to force-free fields with no geometrical symmetries tangent to spheres. 

In the second example ( Section 14. we build exact dynamic isotropic and anisotropic 
plasma equilibrium configurations with magnetic fields tangent to ellipsoids. We start from 
finding a set of non-trivial "vacuum" magnetic fields tangent to ellipsoids, which are then 
transformed into families of dynamic isotropic and anisotropic plasma equilibria by virtue of 
Bogoyavlenskij symmetries ()23|) and " anisotropizing" transformations (|24|) . These solutions 
are well-defined and have a finite magnetic energy in half-space. They model solar coronal 
flares near the active regions of the Sun photosphere. The resulting anisotropic model is an 
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essentially non-symmetric, unlike other available models (see |32j.) It reproduces the features 
of solar flares known from observations, including the presence of thin current sheets. 

In the third example (Section l4.2|) . the coordinate representation is used to build a particu- 
lar trivial "vacuum" magnetic field in prolate spheroidal coordinates. From that field, also by 
the symmetries (|23j) and the transformations (J24j) . we construct families of non-degenerate 
(and generally non-symmetric) isotropic and anisotropic plasma equilibria with dynamics, 
which model the quasi-stationary phase of mass exchange between two spheroidal objects. 
Plasma domains of different geometry and topology may be chosen for the model. 

In the subsection 14.31 conformal coordinate transformations (x, y, z) — > (u, v, w) are em- 
ployed to generate orthogonal coordinates (u,v,w), where exact plasma equilibria can be 
constructed using formulas from Section 13.21 A family of plasma equilibrium configurations 
with non-circular cylindrical magnetic surfaces and realistic values of plasma parameters is 
obtained. 

The suggested approach can be used for the construction of fully 3-dimensional static 
and dynamic plasma equilibrium solutions for MHD and CGL continuum plasma models. 
Configurations with and without geometrical symmetries, for plasma domains of different 
shapes, can be found. 

The author thanks Dr. Kayll Lake and Dr. Oleg Bogoyavlenskij for discussion, and the 
Natural Sciences and Engineering Research Council of Canada (NSERC) for research sup- 
port. 
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